Stability of multidimensional black holes: complete numerical analysis 
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We analyze evolution of gravitational perturbations of D-dimensional Schwarzschild, Reissner- 
Nordstrom, and Reissner-Nordstrom-de Sitter black holes. It is known that the effective potential 
for the scalar type of gravitational perturbations has negative gap near the event horizon. This gap, 
for some values of the parameters Q (charge), A (cosmological constant) and D (number of space- 
time dimensions), cannot be removed by S-deformations. Thereby, there is no proof of (in)stability 
for those cases. In the present paper, by an extensive search of quasinormal modes, both in time and 
frequency domains, we shall show that spherically symmetric static black holes with arbitrary charge 
and positive (de Sitter) lambda-term are stable for D = 5, 6, ... 11. In addition, we give a complete 
numerical data for all three types (scalar, vector and tensor) of gravitational perturbations for multi- 
dimensional black holes with charge and A-term. The influence of charge, A-term and number of 
extra dimensions on black hole quasinormal spectrum is discussed. 

PACS numbers: 04.30.Nk,04.50.+h 



I. INTRODUCTION 



With the development of the string theory and higher- 
dimensional brane- world models, black holes in more 
than four space-time dimensions became a key subject in 
the modern high energy physics. Nevertheless, one of the 
first D-dimensional generalizations of the Schwarzschild 
metric was done, yet, in 1963 by Tangherlini p}. This 
solution has attracted considerable interest recent years, 
in the context of brane- world models [2| , and, for asymp- 
totically anti-de Sitter case, within AdS/CFT correspon- 
dence Q- Different generalizations of the Tangherlini 
metric were found [4(. and their properties were in- 
vestigated in One of the most important property 
of a black hole metric is stability. Space-times unsta- 
ble against small metric perturbations should not ex- 
ist, or need some other, probably quantum, mechanisms 
to provide stability. If a black hole is unstable, small 
perturbations must grow with time, being governed by 
characteristic complex frequencies lo, called quasinor- 
mal modes. These frequencies form a spectrum of the 
response of a black hole to the external perturbation. 
Contribution of all modes can be seen in the time do- 
main [6j. Until now, during recent few years, evolu- 
tion of perturbations of the D-dimensional black holes 
was investigated in a number of works. A lot of pa- 
pers were devoted to perturbations of scalar test fields 
propagating in the vicinity of spherically symmetric black 
holes 7] . Scalar QNMs of higher-dimensional Kerr black 
holes were analyzed in [§j. Quasinormal spectrum of 
Standard Model fields for brane-localized black holes 
were investigated in [§, [n| El, EH- In [Uj a variety 
of brane-localized backgrounds (Schwarzschild, Reissner- 
Nordstrom and Schwarzschild-(anti) de Sitter) were in- 



vestigated, and it was shown that an increase in the 
number of hidden, extra dimensions results in the faster 
damping of all fields living on the brane. In [12] rotat- 
ing black holes on the brane were considered, and, it 
was shown that the brane-localized field perturbations 
are longer-lived when the higher-dimensional black hole 
rotates faster. 

In the context of Tev-gravity scenarios, the gravita- 
tional energy loss in high energy particles collisions was 
considered in [l3[, where was found energy spectra, to- 
tal energy and angular distribution of the emitted grav- 
itational radiation. The further steps were analysis of 
the graviton emission form the higher dimensional black 
holes [14j and the close limit analysis for head-on collision 
of two black holes in higher dimensions [HI, [l(| • 

In addition, some analytical solutions were found in 
different D > 4 space-times [ljj]- The high overtone 
asymptotics for perturbations of D-dimensional black 
holes were considered in [l8[ . 

After appearing of the paper of Kodama and Ishibashi 
(l9l I20I l2l| , where gravitational perturbations of the D- 
dimensional spherically symmetric black hole with charge 
and A-term were reduced to the Schrodinger wave-like 
equations: 



drl 



+ {J 2 - Vi(r))^i = 0, 



(1) 
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the low-laying gravitational quasinormal modes were 
found found in [22J with the help of the WKB method 
[H|]. In (P), i runs T, V, S, because for D > 4, in addi- 
tion to the well-known in four dimensions axial and polar 
perturbations, which are vector and scalar types of grav- 
itational perturbations respectively, there appears a new 
type of gravitational perturbations called tensor one. For 
D-dimensional Reissner- Nordstrom- (anti)-de Sitter black 
holes the tensor type of gravitational perturbations coin- 
cide with perturbations of test scalar filed. In the Ein- 
stein theory with the Gauss-Bonnet term, from works 
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[2J, |25, 26] , it follows that this coincidence does not take 
place. 

An important feature of the gravitational perturba- 
tions is that the effective potentials which governs the 
scalar type of the perturbations has negative gap for 
D > 7, and this gap becomes deeper when D is increased. 
If the effective potential V*(r) is positive definite every- 
where outside the black hole, event horizon, the differen- 
tial operator 

— + UJ 2 
dr 2 

is positive self-adjoint operator in the Hilbert space of 
the square integrable functions of r* . As a result, any 
solution of the wave equation with compact support is 
bounded, what implies stability. In some cases this neg- 
ative gap of the potential can be removed by the so-called 
S-deformation of the potential which does not change 
the property of stabili ty [2l| . In this way Ishibashi 
and Kodama in [l9|, HtJ Hl[ , snowed that D-dimensional 
Schwarzschild black holes are stable against all three 
types of gravitational perturbations. Yet, the potential 
for scalar type of gravitational perturbations of Reissner- 
Nordstrom black holes cannot be made positive definite 
by S-deformation, and therefore the question of stability 
of charged multidimensional black holes left open. That 
is why to judge about (in)stability of black holes we have 
to consider evolution of gravitational perturbations of D- 
dimensional black hole with special emphasis in scalar 
type of gravitational perturbations. 

Note, that with the help of the WKB method, one 
can find only low laying modes for not very large D [13] • 
Therefore we have to apply here the Frobenius method 
[27l ]. Yet, even Frobenius method shows rather bad con- 
vergence for higher D and, this made us to apply, alter- 
natively, the time-domain approach. In the present paper 
we have found all three types of gravitational quasinor- 
mal modes and have shown that the D=5-ll Reissner- 
Nordstrdm(-de Sitter) black holes are stable for any val- 
ues of charge or A-term. 

The paper is organized as follows. In Sec. [Til we 
give the basic formulas and discuss the range of param- 
eters for the D-dimensional Reissner-Nordstr6m(-de Sit- 
ter) metric and of its gravitational perturbations. Sec 
Mil describes the numerical method we used in time and 
frequency domain. Sec IIIII relates the obtained quasi- 
normal modes and time-domain profiles for Reissner- 
Nordstr6m(-de Sitter) black holes in D = 5, 6,... 11 
space-time dimensions. Finally we discuss the obtained 
results and conclude that all considered cases are stable. 



II. BASIC FORMULAS 

The metric of the D-dimensional Reissner-Nordstrom- 
de-Sitter black holes is given by the formula 



dr 2 

ds 2 = f(r)dt 2 - — - - r 2 (dd 2 + sin 2 6d<p 2 ). (2) 
f( r ) 



where 



f(r) = 1 — X + Z — Y, 



(3) 



2M 2Ar 2 
X = ^-, Y = 



rd-l 



z = 



Q 2 



,2d-2 ' 



d(d+l) r* 
The wave equation can be re-written in the form 
d 2 d 2 \ 

W -9^)nt,r) = -V(rm,r), (4) 

where the tortoise coordinate r* is defined as 

dr 



dr* 



The effective potentials for tensor 



V T {r) 



f(r) 



X + d + 



f'(r)rd , f(r)d(d-2) 



2 4 
and vector gravitational perturbations 



(5) 



(6) 



Vv(r) = ^(\ + d±yX 



(d 2 - I) 2 
2 1 ' +2Z\d(d-l) 



d{d - 2)(1 —Y) + Zd(hd - 2) - X(d 2 + 2) 



(7) 



are designated as Vr, Vv respectively. 

The effective potential for the scalar type is given by 



where 



H{r) = A + X 



^+1) 



(8) 



(9) 
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U(r) = (-d 3 (d + l) 2 {d + 2)X 2 + 4d 2 (d + 1) (d(d 2 + 6d - A)Z + 3(d - 2)A) X - 12d 5 (3d - 2)Z 2 - 

-8d 2 (lld 2 - 26d + 12) XZ - 4(d - 2)(d - 4)A 2 ) Y + d 4 (d + 1) 2 X 3 + d(d + 1) (-3d 2 (5d 2 - 5d + 2)Z+ 
+4(2d 2 - 3d + 4)A + d(d - 2)(d - 4)(d + 1)) X 2 + 4d (d 2 (4d 3 + 5d 2 - lOd + 4)Z 2 - 
- (d(34 - 43d + 21d 2 )A + d 2 (d + l)(d 2 - lOd + 12)) Z - 3(d - 4)A 2 - 3d(d + l)(d - 2)A) X - 
-4d 5 (3d - 2)Z 3 + 4)A + 12d 2 (2(-6d + 3d 2 + 4)A + d 2 (3d - 4)(d - 2)) Z 2 + 4(13d - 4)(d - 2)A 2 Z + 
+8d 2 (lld 2 - 18d + 4)XZ + 16A 3 + 4d(d + 2)A 2 , 



Here we used d = D — 2, A = (I + d)(l — 1) and I is the 
multipole number. 

We shall measure all the quantities in the units of the 
horizon radius. Thus we define 



2M = 1 + Q 2 - 2A/d(d+ 1). 



(10) 



In order to parameterize the cosmological constant, we 
introduce the parameter p = l/r c < 1, where r c is the 
cosmological horizon. The value p = corresponds to 
the pure D-dimensional Reissner-Nordstrom black hole. 
Then, A in terms of Q and p can be find from the equation 



/(r c ) = /(1/p) = 0. 
Solving the equation (jllj) . one can find 
d(d+l) 1 - p 1 ^ + Q 2 (l - 



A 



1-P 



-d-1 



(11) 



(12) 



It is clear that for each A there is some extremal charge 
Qext- The value of Q ex t can be found from 



/'(I) = o. 



(13) 



For any charge Q > Q ex t, the definition ([10]) is not cor- 
rect because in this case, the real event horizon has radius 
greater than one. 

Substituting (TT2")) into (TIB")) , one can find 

(d - 1)(1 - p d+1 ) - p 2 {d + 1)(1 - pf*- 1 ) 

^ CXt (d - 1)(1 - /+ 1 ) - p d +!(d + 1)(1 - p d ' r ) ' 

(14) 

Within this paper we shall use p and q = Q/Q ex t as 
the black hole parameters. Note that the formulas (fT4|) 
and (TT2"| allow us to find Q and A in terms of < p < 1 
and < q < 1. In the particular asymptotically flat case 
(A = 0), we have Qe^t = 1 and, therefore, Q = q. 



III. NUMERICAL METHODS FOR 
FREQUENCY AND TIME DOMAINS 

A. Frequency domain 

By definition, the quasi-normal modes are eigenvalues 
of u> with the boundary conditions which correspond to 




FIG. 1: Effective potentials for gravitational perturbations of 
scalar type, D = 5 (blue)...D = 11 (red) {I = 2, Q = 0, 
A = 0). For higher D both the peak and the negative gap of 
the potential increase. 



the outgoing wave at the cosmological horizon (or spatial 
infinity when p = 0) and the ingoing wave at the black 
hole horizon. 

In order to find quasi-normal spectrum for the consid- 
ered types of perturbations, we shall use the continued 
fraction method [l?) • The equation ([1]) has a different 
kind of singularity at r — 1/ p: if p = 0, the singularity is 
irregular at infinity as usually for the Schwarzschild case 



*(r 



oo oc e 



while for p > the point of the cosmological horizon is a 
regular singularity [28| 



* r 



"//'(£) 



It turns out, that we can simply exclude this singularity 
by the introducing of the new function 



where $ satisfies 

d 2 $ 



dr* 



2iuj 



d$ 

dr* 



V4> = 0. 



(15) 




FIG. 2: Effective potentials for gravitational perturbations 
of scalar type, Q = (blue), Q = 0.2 (light blue), Q = 0.4 
(green), Q = 0.6 (light green), Q = 0.8 (yellow), Q = 0.98 
(red) (I = 2, D = 10, A = 0). Increasing of the charge Q 
cause the negative gap to move upwards. For some Q the 
minimum of the potential becomes positive. 




FIG. 3: Effective potentials for gravitational perturbations of 
scalar type, I — 2,3,4,5 (blue, green, yellow, red) (D = 11, 
Q = 0, A = 0). For high multipole numbers the potential 
minimum disappears. Thus for 8 < D < 10 the negative gap 
exists only for the lowest multipole number 1 = 2. 




FIG. 4: Time-domain profiles for gravitational perturbations 
of scalar type (Q = 0, A = 0) for D = 5 (blue). . . D — 11 
(red) at the same point r = 2. Profile for higher D decays 
quicker. 




FIG. 5: Time-domain profiles for gravitational perturbations 
of scalar type (Q = 0, D = 11) for p = 0.3 (blue), p = 0.5 
(green), p = 0.7 (yellow), p = 0.8 (orange), p — 0.9 (red). 
Profile for higher p decays slower. 
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FIG. 6: Time-domain profiles for gravitational perturbations 
of scalar type (A = 0, D = 11) for Q = (dark blue), Q = 
0.5 (blue), Q = 0.6 (light blue), Q = 0.7 (green), Q = 0.8 
(yellow), Q = 0.9 (orange), Q = 0.98 (red). Profile for higher 
Q decays slower. 



In order to satisfy the QNM boundary conditions, the 
function $ has to be convergent at r = 1/ p. Thus, the 
only valuable singularity remains at r = 1. Taking into 
account that there is only the ingoing wave at the event 
horizon, one can find 

^(r-ajcxCr-l)- 2 ^'^. 
Finally, we obtain 

1\ 



$(r) 



1 - - 

r 



y(r), 



(16) 



where y(r) is convergent on r G [1, 1/p]. The appropriate 
Frobenius series for y(r) are 



y{ r ) = J2 art 



r -11- Rp 
r-R 1- p 



(17) 



where R < 1 is an arbitrary parameter. 



FIG. 7: Real and imaginary part of the fundamental quasi- 
normal frequency of gravitational perturbations of scalar type 
as a function of q and p (D = 5, I = 2). 



FIG. 8: Real and imaginary part of the fundamental quasi- 
normal frequency of gravitational perturbations of scalar type 
as a function of q and p (D = 8, I — 2). 
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The expansion (fT7|) is valid if all singular points (except 
those on the event and the cosmological horizons) satisfy 
inequality 



11- Rp 



R 1 



> 1. 



(18) 



Since for the effective potential of the gravitational per- 
turbations of vector and tensor types the singular points 
are 1) r = and 2) the solutions of the equation f(r) = 0, 
we can always find such < R < 1 to satisfy inequality 
(fT8| for all the additional singular points. The latter fol- 
lows from the definition of the horizons, which are the 
two larger positive solutions of the equation f(r) = 0. 

After R is fixed, following [29[ , one can find the recur- 
rence relation for the coefficients a n in (fP7|) and use the 
continued fraction method with the generalized Nollert 
improvement (30j . It turns out, that the convergence of 
the method becomes worse when R — ► 1. Thus one has to 



choose the minimal value of R that satisfies the inequal- 
ity (fT5|) for all singular points, except those at the event 
and the cosmological horizons. For D < 9 this point is 
R = r_, which is the remaining positive root of f(r) (or 
R = if f(r) has no more positive root), and the contin- 
ued fraction method converges without the generalized 
Nollert improvement. For larger D such R is, generally, 
not a singular point of the equation (|15[) and the conver- 
gence of the method is provided by the procedure [29l |. 
Alternatively, for D > 9 one can continue the Frobenius 
series through some midpoints. This method was used in 
[3l| and our results are completely in agreement with it. 

For the scalar- type gravitational perturbations, the 
equation (|15[) has additional singularities, which are the 
roots of H(r). Thus, in higher dimensions D > 7, for 
some values of p, q and I there is no such R < 1, which 
satisfies the inequality (fl8|) for all the singular points. 
Therefore, for those cases one should use only the tech- 
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FIG. 9: Real and imaginary part of the fundamental quasi- 
normal frequency of gravitational perturbations of scalar type 
as a function of I (D = 9, q = 0.5, p = 0.3). The solid line 
represents high multipole WKB formula. 
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nique developed in [3l| . Practically, the continued frac- 
tion coefficients appear to be so complicated, that we are 
unable to compute QNMs with this method during rea- 
sonable time, even for low dimensions, except the case 
when (0 = 0. 



B. Time domain 

Alternatively, we shall use the time-domain analysis 
based on the standard integration scheme for the wave- 
like equations ([J|, which is described for instance in Q. 
This approach has no limitations of the parameters for 
the scalar-type perturbations. Moreover, it appears to 
be quicker for the perturbations of scalar type, when D 
is large or p ^ 0. Since we can see the real signal which 
contains all the frequencies, this method also provides 
the most direct evidence of the black hole stability. It 
allows also to find out asymptotical tails. 

In detail, we applied a numerical characteristic integra- 
tion scheme, that uses the light-cone variables u = t — r * 
and v = t+r+. In the characteristic initial value problem, 
initial data are specified on the two null surfaces u = uq 



and v — vq. The discretization scheme we used, is 

_ a ,v(wMw) + v(eme) +0{a4)Ai9) 

8 

where we have used the following definitions for the 
points: N = {u+A, v+A), W = (u+A, v),E= (u, v+A) 
and S = (it, v). 

IV. EVOLUTION OF GRAVITATIONAL 
PERTURBATIONS 

A. Scalar type of gravitational perturbations 

It is essential for our computations that the effective 
potential for charged black holes have the negative gap 
which cannot be "removed" by S-deformation [l^. I20L l2lj| 
for D > 8. This gives us the cases for suspecting instabil- 
ity. The fundamental quasinormal modes for scalar type 
of gravitational perturbations computed with Frobenius 
and time-domain techniques are given in Tables [IJ HT1 
For D = 5, 6, 7 Reissner- Nordstrom black holes (A = 0), 
the real parts of QNMs increase with charge until some 
maximum, then decrease. For D > 8 Reu, for small Q, 
decreases as a function of Q, then increases and reaches 
some maximum and then, again decreases for close to the 
extremal Q. The imaginary part of the frequencies as a 
function of charge Q monotonically decreases for all D 
(Table HI). Note that in table U the calculations of QNMs 
were done for D = 5, 6, 7 with both time-domain and 
Frobenius approaches. Yet, for higher D, the Frobenius 
method needs considerable computer time for converging. 
The time-domain picture of the decay for black holes with 
high D = 10, 11, . . . are different from lower dimensional 
cases: if for D — 5, 6, ... 9 we see always regular picture 
of time-domain decay (Fig. 4), for D — 10, 11, and for 
small Q and A the ringing is more irregular (in fact, the 
amplitude is irregular, that probably means the superpo- 
sition of a few dominating modes) and quickly changes 
into tail behavior (Fig. 5, Fig. 6). The time domain pic- 
ture becomes again regular at larger Q and A (Fig. 5, Fig. 
6). 

The fundamental modes for Schwarzschild -de Sitter 
black holes (Q = 0) are given in Table [Til The real os- 
cillation frequency usually monotonically decreases as a 
function of A-term. Note that literature on QNMs of the 
scala r type o f gravitational perturbations in D > 4 is very 
poor [22|,|32[. Our computations show perfect coincidence 
with Frobenius method data of I32H and very good coinci- 
dence with the WKB results of [22| for scalar type of grav- 
itational perturbations. Alternative approach for find- 
ing gravitational quasinormal modes is to consider the 
head-on collisions for higher dimensional theories (l5l.[l^|. 
The comparison with our accurate numerical values of 
QNMs gives very good agreement for not large D: Thus 
for D = 5 and 6, we have w = 0.94773 - 0.25609i, 
lo = 1.13690 — 0.30357i respectively, while paper [l5[ give 
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TABLE I: Fundamental modes (I = 2, A = 0) of gravitational perturbations of scalar type calculated by the continued fraction 
method (the first value) and by the time-domain simulation (the second value). For higher D the behavior of the time-domain 
profile is very irregular. Thus QNMs finding with high accuracy is more complicated. 
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TABLE II: Fundamental modes (I = 2, Q = 0) of gravitational perturbations of scalar type calculated by the time-domain 
simulation. For higher D the behavior of the time-domain profile is very irregular. Thus QNMs finding with high accuracy is 
more complicated. 
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1.31868 - 


0.39970* 


1.53691 - 





59412* 


1.95602 - 


0.86639* 


2.41 


- 0.99* 


0.2 


0.89204 - 


0.25366* 


1.07169 - 


0.30664* 


1.25828 - 


0.39567* 


1.45301 - 





56754* 


1.81939 - 


0.87158* 


2.28 


- 1.00* 


0.3 


0.82569 - 


0.24749* 


0.99203 - 


0.30681* 


1.16108 - 


0.38781* 


1.33119 - 





52832* 


1.56112 - 


0.84743* 


2.04 


- 1.01* 


0.4 


0.73774 - 


0.23477* 


0.88551 - 


0.30073* 


1.03263 - 


0.37458* 


1.17909 - 





48504* 


1.31893 - 


0.67923* 


1.70 


- 1.01* 


0.5 


0.63289 - 


0.21293* 


0.75787 - 


0.28372* 


0.87872 - 


0.35286* 


0.99978 - 





44032* 


1.11688 - 


0.56723* 


1.209 


-0.801* 


0.6 


0.51662 - 


0.18156* 


0.61591 - 


0.25212* 


0.70774 - 


0.31741* 


0.79938 - 





38928* 


0.89080 - 


0.47868* 


0.978 


- 0.60H 


0.7 


0.39236 - 


0.14257* 


0.46597 - 


0.20493* 


0.52888 - 


0.26326* 


0.58941 - 





32283* 


0.64959 - 


0.38903* 


0.710 


- 0.468* 


0.8 


0.26341 - 


0.09797* 


0.31181 - 


0.14428* 


0.34978 - 


0.18904* 


0.38262 - 





23377* 


0.41294 - 


0.28050* 


0.442 


- 0.332* 


0.9 


0.12993 - 


0.04982* 


0.15447 - 


0.07427* 


0.17045 - 


0.09844* 


0.18471 - 





12283* 


0.19478 - 


0.14745* 


0.203 


-0.173* 



the values w = 0.947 - 0.256*, to = 1.139 - 0.305*, i.e. 
the error in [l5| is indeed very small and is less than 1 
percent. For D = 9, 10, 11 the method in [l5[ shows su- 
perposition of the two dominating modes, one of which is 
very close to that obtained here. For instance for D = 4, 
they have u> = 2.47 — 0.99*, while we obtained with time- 
domain approach li — 2 Ah — 0.98*. This small differ- 
ence between our results is quite natural because differ- 
ent approximations were used: we considered the regime 
of small perturbations and in [loT ] a close limit approxi- 



mation was used. Therefore we conclude that the close 
limit approximation with Brill-Lindquist initial data is 
indeed a very accurate approach. 

When both parameters Q and A are not vanishing, the 
quasinormal frequencies are again damping. An example 
of behavior of real and imaginary of w as a function of q 
and p (which parameterize all the values of Q and A) is 
given in Fig [7] for D = 5 case and Fig [8] for D = 8 case. 

At large multi-pole numbers £, the WKB approach [23| 
becomes very accurate and one can obtain usually an 
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analytical formula for the QN frequencies. For our case of 
Reissner-Nordstrom-de Sitter black holes, a cumbersome 
form of the effective potential with a lot of parameters 
makes it impossible to find a short QNM formula, which 
could be easily interpreted for general values of M, Q, 
A, D. Yet, if we fix all the above parameters the WKB 
formula in the lowest order gives, for example for D = 9, 
Q = 0-5Q ex u P = 0.3 (see Fig El): 

w = 0.60897(^ + 3) - 1.42338(n+ l/2)i, £ -> oo (20) 

Here we can see that in a similar fashion with a four- 
dimensional case, the real oscillation frequency is pro- 
portional to £ and the damping rate is approaching some 
constant value. Therefore one should not expect any in- 
stability from the regime of high £. Indeed, higher £ just 
raise the potential barrier and serve as a stabilizing fac- 
tor: perturbations with higher £ are more stable than 
perturbations with lower £. 

Finally, summarizing the extensive search for quasinor- 
mal modes for scalar type of gravitational perturbations, 
we conclude that the response to the external perturba- 
tions of the Reissner-Nordstrom-de Sitter black holes are 
represented by damping oscillations and therefore for all 
values of charge and A-term in D = 5,6, ...11 space- 
time dimensions, the Reissner-Nordstrom-de Sitter black 
holes are stable. 



B. Vector and tensor types of gravitational 
perturbations 

The existing literature on vector and tensor type 
of gravitational perturbations is much fuller for the 
Schwarzschild case [22|, [3l|, [HJ , . In addition the effective 
potential is positive definite for these types of pertur- 
bations, so we do not suspect instability in them. The 
fundamental modes for vector type " + " gravitational 
perturbations are given in Table IIIII the real oscilla- 
tion frequency of the dominating mode is growing as the 
charge Q grows, and after reaching some maximum, is 
decaying. For D = 5, the maximum of Reu happens 
for Q » Q.bQ ex tr- This value of charge, at which Reuj 
has maximum, is increasing for higher D, and is about 
Q « 0.7Q ex tr for D = 10. The imaginary part of the 
frequency is monotonically decreasing as Q is increasing 
for vector type " + " perturbations for all D. 

The vector type " — " of gravitational perturbations 
for Reissner-Nordstrom case has rather different behav- 
ior (see Table IIII[) : both Reto and Imuj are monotoni- 
cally decaying as a function of Q. Let us note that for 
vector type of gravitational perturbations we perfectly 
reproduce the results of the pure Schwarzschild limit of 
the previous papers [H, [3l|, H2] . Thus, for example, for 
D = 5we get in Table Ell w = 1.13400-0.32752i for vec- 
tor type " — " of gravitational perturbations obtained with 
the time-domain approach, while the Frobenius method 
in [32j give u> = 1.1340 — 0.3275i. Let us remember that 



for vector type of gravitational perturbations there are 
the two types " — " and " +" of perturbations, because in 
fact we have the perturbations of the Maxwell field and of 
gravitational field, which are coupled. Thereby the type 
" + " of vector perturbations reduces to the perturba- 
tions of the test Maxwell field when the electromagnetic 
background vanishes Q = 0. 

Dependence of the QNMs on the values of the A-term 
parameter is the same for both types of vector perturba- 
tions: both Imuj and Rem are monotonically decreasing 
when A is increasing. 

The tensor type of gravitational perturbations is 
known to coincide with perturbations of test scalar field 
fl9l |. and the corresponding quasinormal spectrum is 
isospectral as well [22j . The fundamental modes for pure 
Reissner-Nordstrom and Schwarzschild-de Sitter cases 
are given in Tables VII and VIII. Both real and imag- 
inary parts of quasinormal frequencies are monotonically 
decreasing as charge or A-term is increasing. 

This way one has quasinormal modes for all types 
of gravitational perturbations analyzed for Reissner- 
Nordstrom-de Sitter black holes. 



V. CONCLUSIONS 

Multidimensional black holes are subject of intensive 
research in brane- world theories during recent years, be- 
cause they are expected to be produced in particles col- 
lisions in Large Hadron Collider or in cosmic showers ex- 
periments. The basic feature one needs to provide relia- 
bility of the black hole metric is stability. Yet the stability 
of charged multidimensional black holes with charge and 
A-term for D > 8 was an open question until the present 
research. We have performed here an extensive time and 
frequency domains investigation of a response of multidi- 
mensional black holes with arbitrary charge and A-term 
and showed that the response is decaying and thereby the 
black holes under consideration are stable. As a by prod- 
uct we presented a detailed data on quasinormal frequen- 
cies for the considered black holes and the dependence of 
the frequencies on such parameters as black hole charge 
and the value of the cosmological term. 

The are a few important points beyond our considera- 
tion: first of all, it is Gauss-Bonnet generalizations of the 
Einstein solutions [25|, [26| , where an instability can hap- 
pen even in asymptotically flat case 33] in D — 5 case. 
There numerical analysis of the gravitational perturba- 
tions of the Gauss-Bonnet black holes for different D is in 
progress [341 ] . Another important and interesting case is 
asymptotically anti-de Sitter charged black holes, when 
the instability also exists [351 ] and can be interpreted as a 
phase transition in the dual temperature conformal field 
theory (36|. 

Apparently, an interesting alternative model for higher 
dimensional black holes as a window in extra dimensions 
is the Kaluza-Klein black holes with squashed horizons 
[37l ]. The gravitational perturbations for these solution 
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TABLE III: Fundamental modes (I = 2, A = 0) of gravitational perturbations of vector("+") type calculated by the time-domain 
simulation. 



Q 




D = 


= 5 




D = 


= 6 




D = 


= 7 




D = 


= 8 




D = 


= 9 




D 


= 10 


n 

u 




I 


46854 — 


n 
u 


OOZ4:Z£ 


I 


97832 — 


n 

u 




2 


47008 — 


u 


fi9fizirw 

UiU'iUt 


2 




n 




3 


43062 — 


n 
u 


8^1 An 


3 


Q055Q 


n 

— u 




n 

u 


1 

± 


1 


47080 — 


n 

u 




1 


98180 — 


n 

u 




2 


47403 — 


n 

u 




2 


95761 — 


u 


74970? 


3 


43537 _ 


u 




3 


91135 


— u 




n 
u 


O 
z 




47742 — 


n 
u 


1A P^l i 


1 


99114 _ 


n 
u 




2 


48556 — 


n 
u 


OlOoo? 


2 


97099 — 


n 
u 


7^/1 Q1 i 

t o^ry i? 


3 


45092 — 


n 
u 




3 


92777 


— u 




n 

u 


Q 

o 


x 


48547 — 


n 

u 


QQRIS79 
OOUU I i 


2 


00371 — 


n 
u 




2 


50181 — 


n 
u 




2 


99037 — 


n 
u 


791 77i 


3 


47300 — 


n 
u 




3 


Q5998 


_ n 


you i ut 





4 


1 


49155 — 





32391? 


2 


01585 — 





46226? 


2 


51897 — 





58770? 


3 


01180 — 





70309i 


3 


49gl9 _ 





81005i 


3 


98080 


-0 


90983i 





5 


1 


49210 - 





30723? 


2 


02337 - 





44241?' 


2 


53244 - 





56546? 


3 


03033 - 





67894? 


3 


52125 - 





78433i 


4 


00792 


-0 


88286i 





6 


1 


48353 - 





28708? 


2 


02189 - 





41860? 


2 


53713 - 





53882? 


3 


04036 - 





65008i 


3 


53610 - 





75362i 


4 


02713 


-0 


85065i 





7 


1 


46249 - 





26477? 


2 


00727 - 





39257? 


2 


52816 - 





50988?' 


3 


03633 - 





61882? 


3 


53655 - 





72040i 


4 


03167 


-0 


81580i 





8 


1 


42669 - 





24279? 


1 


97680 - 





36751?' 


2 


50236 - 





48241? 


3 


01453 - 





58938? 


3 


51835 - 





68929i 


4 


01676 


-0 


78326i 





9 


1 


37633 - 





22431?' 


1 


93086 - 





34705? 


2 


46010 - 





46053? 


2 


97522 - 





56631? 


3 


48165 - 





66523i 


3 


98235 


-0 


75837i 





98 


1 


32775 - 





21336? 


1 


88525 - 





33525?' 


2 


41682 - 





44832?' 


2 


93363 - 





55377?' 


3 


44148 - 





65248i 


3 


94338 


-0 


74548i 



TABLE IV: Fundamental modes (I — 2, A = 0) of gravitational perturbations of vector(" — ") type calculated by the time-domain 
simulation. 



Q 




D = 


= 5 




D = 


= ( 






D = 


= 7 




D = 


= 8 




D = 


= 9 




D 


= 10 







1 


13400 - 





32752i 


1 


52467 - 





474 12?' 


1 


93447 - 





61230?' 


2 


35885 - 





73871?' 


2 


79285 - 





85412? 


3 


23345 


-0 


96025? 





1 


1 


12873 - 





32542? 


1 


51974 - 





47195?' 


1 


93008 - 





60995? 


2 


35455 - 





73593? 


2 


78882 - 





85126? 


3 


22923 


-0 


95699? 





2 


1 


11323 - 





31965? 


1 


50551 - 





46553? 


1 


91694 - 





60268? 


2 


34200 - 





72777? 


2 


77645 - 





84240? 


3 


21691 


-0 


94756* 





3 


1 


08935 - 





31042i 


1 


48324 - 





45520? 


1 


89595 - 





59098?' 


2 


32163 - 





71475?' 


2 


75614 - 





82835? 


3 


19652 


-0 


93271? 





4 


1 


05895 - 





29829i 


1 


45430 - 





44154?' 


1 


86800 - 





57559?' 


2 


29396 - 





69786? 


2 


72819 - 





81033? 


3 


16820 


-0 


91389? 





5 


1 


02349 - 





28402? 


1 


41985 - 





42551?' 


1 


83391 - 





55781?' 


2 


25963 - 





67875i 


2 


69314 - 





79031? 


3 


13248 


-0 


89327? 





6 





98411 - 





26862? 


1 


38094 - 





40854? 


1 


79473 - 





53951? 


2 


21974 - 





65967? 


2 


65228 - 





77079? 


3 


09083 


-0 


87359? 





7 





94199 - 





25339? 


1 


33903 - 





39228?' 


1 


75230 - 





52260?' 


2 


17659 - 





64270? 


2 


60824 - 





75390? 


3 


04653 


-0 


85697? 





8 





89871 - 





23937? 


1 


29619 - 





37763?' 


1 


70919 - 





50774?' 


2 


13312 - 





62808?' 


2 


56415 - 





73951? 


3 


00213 


-0 


84311? 





9 





85578 - 





22668? 


1 


25410 - 





36428? 


1 


66671 - 





49423? 


2 


09047 - 





61471?' 


2 


52102 - 





72627? 


2 


95869 


-0 


83045? 





98 





82248 - 





21747? 


1 


22123 - 





35412? 


1 


63310 - 





48428? 


2 


05718 - 





60463?' 


2 


48645 - 





71566? 


2 


92469 


-0 


82086?' 



would be especially actual, as they would probably show 
gravitational instability of the Gregory-Laflamme type. 
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TABLE V: Fundamental modes (I = 2, Q — 0) of gravitational perturbations of vector(" + ") type calculated by the continued 
fraction method (the first value) and by the time-domain simulation (the second value) . 



p 




D = 




) 




D = 


= 6 




D = 


= 7 




D = 


= 8 




D = 


= 9 




D 


= 10 







1 


46854 - 





35243i 


1 


97830 - 





49642i 


2 


47006 - 





62642i 


2 


95278 - 





745324 


3 


43056 - 





855184 


3 


90550 


-0 


957494 






1 


46854 - 





35242i 


1 


97832 - 





496414 


2 


47008 - 





626404 


2 


95282 - 





745294 


3 


43062 - 





855144 


3 


90559 


-0 


957434 





1 


1 


44629 - 





347524 


1 


95277 - 





491234 


2 


44074 - 





621134 


2 


91936 - 





740054 


3 


39283 - 





849984 


3 


86332 


- 


952414 






1 


44629 - 





3475 li 


1 


95278 - 





491224 


2 


44076 - 





621124 


2 


91940 - 





740024 


3 


39289 - 





849944 


3 


86341 


-0 


952354 





2 


1 


38107 - 





3329H 


1 


87596 - 





475 144' 


2 


35219 - 





604494 


2 


81847 - 





723264 


3 


27906 - 





833274 


3 


73627 


-0 


935914 






1 


38107 - 





3329H 


1 


87597 - 





475 1 34' 


2 


35221 - 





60447i 


2 


81850 - 





723234 


3 


27911 - 





833234 


3 


73635 


-0 


935864 





3 


1 


27723 - 





30 9 084' 


1 


74877 - 





447024 


2 


20360 - 





574334 


2 


64868 - 





692144 


3 


08771 - 





801754 


3 


52289 


-0 


904334 






1 


27724 - 





30 9 084' 


1 


74877 - 





4470H 


2 


20361 - 





574314 


2 


64871 - 





692124 


3 


08776 - 





801724 


3 


52296 


-0 


904294 





4 


1 


14122 - 





276984 


1 


57512 - 





406124 


1 


99630 - 





528094 


2 


40957 - 





642734 


2 


81743 - 





750444 


3 


22147 


- 


851884 






1 


14122 - 





276984 


1 


57513 - 





406114' 


1 


99631 - 





528084 


2 


40959 - 





642714 


2 


81747 - 





750424 


3 


22152 


-0 


851854 





5 





98025 - 





2382H 


1 


36231 - 





353134 


1 


73596 - 





46457i 


2 


10458 - 





571704 


2 


46965 - 





674134 


2 


83196 


-0 


771824 









98025 - 





23821i 


1 


36231 - 





353134 


1 


73597 - 





464564 


2 


10460 - 





571694 


2 


46967 - 





674124 


2 


83199 


- 


771794 





6 





80110 - 





194664 


1 


11935 - 





290564 


1 


43247 - 





385474 


1 


74298 - 





478894 


2 


05196 - 





57025i 


2 


35987 


-0 


659124 









80110 - 





19466i 


1 


11935 - 





290554 


1 


43248 - 





38547i 


1 


74299 - 





478884 


2 


05198 - 





570244 


2 


35989 


-0 


659104 





7 





60947 - 





14800i 





85494 - 





221644 


1 


09744 - 





295264 


1 


33861 - 





368814 


1 


57929 - 





442114 


1 


81992 


-0 


51489i 









60946 - 





148014 





85494 - 





221644 


1 


09744 - 





295264 


1 


33862 - 





368814 


1 


57930 - 





442104 


1 


81993 


-0 


514884 





8 





40990 - 





099474 





57648 - 





149154' 





74150 - 





198904 





90581 - 





248804 


1 


06987 - 





298874 


1 


23394 


-0 


349094 









40990 - 





099474 





57648 - 





14915i 





74150 - 





198904 





90581 - 





248804 


1 


06987 - 





298874 


1 


23395 


-0 


349094 





9 





20588 - 





049944 





28997 - 





074914' 





37338 - 





099894 





45650 - 





124904 





53948 - 





149954 





62239 


-0 


175054 









20588 - 





049944' 





28997 - 





074914' 





37338 - 





09989i 





45650 - 





124904 





53948 - 





149954 





62239 


-0 


175054 



TABLE VI: Fundamental modes (I = 2, Q = 0) of gravitational perturbations of vector(" — ") type calculated by the continued 
fraction method (the first value) and by the time-domain simulation (the second value) . 



p 




D = 


_ 1 






D = 


= 6 




D = 


= 7 




D = 


= i 


i 




D = 


= 9 




D 


= 10 







1 


13400 - 





32 7524' 


1 


52466 - 





474124 


1 


93446 - 





612314 


2 


35883 - 





738744' 


2 


79281 - 





854154 


3 


23339 


- 


960314 






1 


13400 - 





32 7524' 


1 


52467 - 





47412i 


1 


93447 - 





612304 


2 


35885 - 





7387H 


2 


79285 - 





854124 


3 


23345 


-0 


960254 





1 


1 


11435 - 





324454 


1 


49999 - 





471204' 


1 


90445 - 





609804 


2 


32356 - 





736704 


2 


75235 - 





852554 


3 


18778 


-0 


959074 






1 


11435 - 





324454 


1 


50000 - 





471194' 


1 


90447 - 





609794 


2 


32358 - 





736684 


2 


75238 - 





852524 


3 


18784 


-0 


959024 





2 


1 


05738 - 





314484 


1 


42675 - 





460854 


1 


81478 - 





600354 


2 


21794 - 





728544 


2 


63114 - 





845614 


3 


05117 


-0 


953214 






1 


05738 - 





314484' 


1 


42676 - 





460844 


1 


81479 - 





600344 


2 


21796 - 





728524 


2 


63117 - 





845584 


3 


05122 


-0 


953164 





3 





96883 - 





296084 


1 


30898 - 





439234 


1 


66798 - 





578634 


2 


04364 - 





708274 


2 


43048 - 





827084 


2 


82473 


- 


936284 









96883 - 





296084 


1 


30898 - 





439234 


1 


66798 - 





578624 


2 


04365 - 





708264 


2 


43051 - 





827064 


2 


51248 


-0 


898124 





4 





85670 - 





268724 


1 


15581 - 





402984 


1 


47213 - 





538114 


1 


80689 - 





667324 


2 


15518 - 





787344 


2 


51245 


-0 


898154 









85671 - 





268724 


1 


15581 - 





402974 


1 


47213 - 





538104 


1 


80690 - 





667314' 


2 


15520 - 





787324 


2 


12471 


-0 


826484 





5 





72870 - 





233524 





97927 - 





352264 


1 


24192 - 





475394 


1 


52193 - 





598194 


1 


81754 - 





715924 


2 


12469 


-0 


826504 









72870 - 





233524 





97927 - 





352254 


1 


24192 - 





475394 


1 


52194 - 





598184 


1 


81755 - 





715914 


1 


68348 


-0 


710834 





6 





59065 - 





192374 





78984 - 





290494 





99434 - 





393414 


1 


21066 - 





499894 


1 


44075 - 





606714 


1 


68347 


-0 


710854 









59065 - 





192374 





78984 - 





290494 





99434 - 





393414 


1 


21066 - 





499884 


1 


44075 - 





606714 


2 


82477 


-0 


936244 





7 





44658 - 





147104' 





59441 - 





22 1 744' 





74231 - 





299394 





89505 - 





380634 


1 


05542 - 





464924 


1 


22478 


-0 


550814 









44658 - 





147104' 





59441 - 





22 1 744' 





74231 - 





299384 





89505 - 





380634 


1 


05543 - 





464924 


1 


22479 


-0 


550814 





8 





29911 - 





099234 





39673 - 





149204 





49222 - 





200254 





58809 - 





252894 





68581 - 





307434 





78648 


-0 


363964 









29911 - 





099234 





39673 - 





149204 





49222 - 





200254 





58809 - 





252894 





68581 - 





307434 





78648 


-0 


363964 





9 





14990 - 





049914 





19842 - 





074914 





24523 - 





100064 





29132 - 





125424 





33718 - 





151104 





38309 


-0 


17716i 









14990 - 





049914' 





19842 - 





074914' 





24523 - 





IOOO64' 





29132 - 





125424 





33718 - 





151104 





38309 


- 


177164 
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TABLE VII: Fundamental modes (I = 2, A = 0) of gravitational perturbations of tensor type calculated by the continued 
fraction method. 



o 


D = 


= 5 


D = 


= 6 


D = 


= 7 


D = 


= 8 


D = 


= 9 


D 


= 10 


o 


1.51057 — 


357547 


9 ni i 53 — 


501 947 

yj . W w -L E?*±fc 


9 49678 — 


631 887 

w. UO -LOO L 


9 97469 — 


750567 
yj • 1 juuut 


3 44889 — 


8601 37 


3 99095 


— 9691 9? 


0.1 


1 50687 — 

-L . w w wo i 


355277 

w . 'J w u i i> 


2 00760 — 


499267 


2 49270 — 


628937 


9 97048 — 


747437 


3.44451 — 


85686? 


3 91654 

O • O -L WWt: 


— 95876? 

W. JUL) 1 W () 


0.2 


1.49575 — 


34858? 


1 qqe;7Q _ 

± .C/tjO 1 is 


491 367 

yj .'rtiJ -L Jut 


9 48041 — 


620307 


9 95782 — 


738997 


3 431 53 — 


84736? 

U . 04: I OU (j 


3 90398 


W. flO^Ot 


0.3 


1.47711 — 


33773i 


1.97601 — 


47872? 

W ■ j: I W 1 


2 45985 — 


60659? 


2 93664 — 


72389? 


3 40980 — 


83247z 


3 88107 


— 93374? 

W • W W J ^tt 


0.4 


1.45084 - 


0.32326i 


1.94822 - 


0.46216? 


2.43101 - 


0.58890?' 


2.90695 - 


0.70552i 


3.37938 - 


0.81368i 


3.85000 


- 0.91466i 


0.5 


1.41684 - 


0.30600?: 


1.91251 - 


0.44296? 


2.39413 - 


0.56881? 


2.86911 - 


0.68503i 


3.34072 - 


0.79303i 


3.81063 


- 0.89397i 


0.6 


1.37526 - 


0.28717?: 


1.86940 - 


0.42278? 


2.35002 - 


0.54830i 


2.82423 - 


0.6646H 


3.29520 - 


0.77285i 


3.76452 


- 0.87411? 


0.7 


1.32676 - 


0.26838i 


1.82018 - 


0.40345? 


2.30048 - 


0.52922?' 


2.77447 - 


0.64603i 


3.24525 - 


0.75485i 


3.71440 


- 0.85670i 


0.8 


1.27290 - 


0.25118i 


1.76696 - 


0.38619? 


2.24787 - 


0.51240? 


2.72233 - 


0.62977i 


3.19345 - 


0.73917i 


3.66296 


- 0.84114? 


0.9 


1.21589 - 


0.23645?: 


1.71184 - 


0.37127i 


2.19317 - 


0.49714? 


2.66790 - 


0.61582i 


3.14034 - 


0.72614z 


3.59458 


- 0.78869? 



TABLE VIII: Fundamental modes (I = 2, Q = 0) of gravitational perturbations of tensor type calculated by the continued 
fraction method. 



p 


D = 


= 5 


D = 


= 6 


D = 


= 7 


D = 


= 8 


D = 


= 9 


D 


= 10 





1.51057 - 


0.35754?: 


2.01153 - 


0.50194?: 


2.49678 - 


0.63188?: 


2.97469 - 


0.75056? 


3.44882 - 


0.86013? 


3.92095 


- 0.96212? 


0.1 


1.48334 - 


0.35388i 


1.98035 - 


0.49836i 


2.46129 - 


0.62843?: 


2.93467 - 


0.74725? 


3.40415 - 


0.85698? 


3.87154 


- 0.95914? 


0.2 


1.40454 - 


0.34210i 


1.88747 - 
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0.61657? 


2.81464 - 


0.73573? 


3.27015 - 
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3.72335 
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0.3 


1.28232 - 


0.32055i 
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2.17943 - 
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2.61539 - 
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3.04733 - 
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0.4 


1.12818 - 
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0.42371i 
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2.73842 - 
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3.13409 


- 0.88045? 


0.5 


0.95396 - 
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0.36922?: 
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0.6 
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0.20040?: 


1.05127 - 


0.30195i 
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0.40307? 


1.61764 - 


0.50244? 
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0.59910? 


2.19479 


- 0.69251? 


0.7 


0.57901 - 


0.15064i 


0.79008 - 


0.22746i 


0.99988 - 


0.30527? 
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0.38368? 
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0.46209? 
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